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Abstract 



We have systematically investigated the semi-inclusive B decays B — * MX, which are manifes- 
tations of the quark decay b — ► Mq, within the framework of QCD-improved factorization. These 
decays are theoretically clean and have distinctive experimental signatures. We focus on a class of 
these that do not require any form factor information and therefore may be especially suitable for 
extracting information on the angles a and 7 of the unitarity triangle. The nonfactorizable effects, 
such as vertex-type and penguin-type corrections to the two-body b decay, b — ► Mq, and hard spec- 
tator corrections to the 3-body decay B — » Mq\q\ are calculable in the heavy quark limit. QCD 
factorization is applicable when the emitted meson is a light meson or a charmonium. We discuss the 
issue of the CPT constraint on partial rate asymmetries. The strong phase coming from final-state 
rescattering due to hard gluon exchange between the final states can induce large rate asymme- 
tries for tree-dominated color-suppressed modes (tt , p°, u)X§. The nonfactorizable hard spectator 
interactions in the 3-body decay B — > Mq\qi, though phase-space suppressed, are extremely impor- 
tant for the tree-dominated modes (tt°, p , ui)Xg, <f>X, J/ipX s , J/ipX and the penguin-dominated 
mode ojX S g. In fact, they are dominated by the hard spectator corrections. This is because the 
relevant hard spectator interaction is color allowed, whereas the two-body semi- inclusive decays 
for these modes are color-suppressed. Our result for B(B — * J/ifjX s ) is in agreement with exper- 
iment. The semi-inclusive decay modes: B° s — > (tt°, p°, uj)X s , p°X sS , — ► {K~X, K*~X) and 
B~ — > (K°X S , K*°X S ) are the most promising ones in searching for direct CP violation. In fact, 
they have branching ratios of order 10 -6 — 10 -4 and CP rate asymmetries of order (10 — 40)%. 
The decays B s — > (ir°, p°, uj)X sS and B~ — > <f>X are electroweak-penguin dominated. Some of them 
have sizable branching ratios and observable CP asymmetries. 
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I. INTRODUCTION 



Recently Beneke, Buchalla, Neubert and Sachrajda (BBNS) JTJ have proposed a QCD- 
improved factorization approach to a class of exclusive S-decays that appears quite promis- 
ing. BBNS suggest that nonfactorizable effects in an exclusive decay B — > M\M 2 with 
recoiled Mi and emitted light meson M 2 are calculable since only hard interactions between 
the (BM\) system and M 2 survive in the heavy quark limit. In this paper we extend the 
application of BBNS idea of QCD factorization to a certain class of semi-inclusive decays. In 
this regard our approach complements the recent works of He et al fHH . The semi-inclusive 
decays that are of special interest originate from the quark level decay, b — > M + q; they 
are theoretically cleaner compared to exclusive decays and have distinctive experimental 
signatures . Since these semi-inclusive decays also tend to have appreciably larger branch- 
ing ratios compared to their exclusive counterparts, they may therefore be better suited for 
extracting CKM-angles and for testing the Standard Model (SM). 

Earlier studies of semi-inclusive decays are based on naive factorization Q or generalized 
factorization H in which nonfactorizable effects are treated in a phenomenological way by 
assuming that, for example, the number of colors is a free parameter to be fitted to the 
data. Apart from the unknown nonfactorizable corrections, this approach encounters another 
major theoretical uncertainty, namely the gluon's virtuality k 2 in the penguin diagram is 
basically unknown, rendering the predictions of CP asymmetries not trustworthy. 

The aforementioned difficulties with the conventional methods can be circumvented in 
the BBNS approach of QCD-improved factorization (QCDF). Indeed, by placing an energy 
cut Em > 2.1 GeV on the meson in the semi-inclusive signal B — > M + X, not only we 
can enhance the presence of the two body quark decays, b — > M + q, but also M then can 
play the role of M 2 and X that of the recoiled meson Mi in the above-mentioned exclusive 
decay B — ► Mi + M 2 , in so far as considerations of BBNS go. Furthermore, a very important 
theoretical simplification occurs in the semi-inclusive decays over the exclusive decays if we 
focus on final states such that M does not contain the spectator quark of the decaying B(B S ) 
meson as then we completely by-pass the need for the transition form factor for B(B S ) — > M. 
Recall that for the exclusive case, in general, we need a knowledge of two such form factors 
if M is a pseudoscalar meson or of four form factors if M is a vector meson. 

The consideration of these semi-inclusive B decays has several other theoretical advan- 
tages over the exclusive ones as well. For one thing, there is no troublesome infrared divergent 
problem occurred at endpoints when working in QCD factorization. As for CP violation, 
contrary to the exclusive hadronic decays, it is not plagued by the unknown soft phases. Con- 
sequently, the predictions of the branching ratios and partial rate asymmetries for B — > MX 
are considerably clean and reliable. 

Recently QCD factorization has been applied to charmless semi-inclusive decays B — > 
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K(K*)X and B — > 0X S in [@,[3|]. In the present paper we will systematically study all semi- 
inclusive decays B -> PX(X S ,X C ), 5 -> VX(X S ,X C ) and S -> J/t/X with P (V) being 
a light pseudoscalar (vector) meson, X the final state containing no charmed or strange 
particles and X s the final state containing a strange quark and X c the final state containing a 
charm quark. We differ from |2],|3| in two main aspects: First, we have included the complete 
twist-3 corrections to the penguin coefficient to be introduced below and electroweak 
penguin-type corrections to the coefficients 07-10 induced by tree 4-quark operators; both 
have been neglected in |2|||. Second, the hard spectator interaction in the 3-body decay 
B — > Mq±q 2 is neglected in [@,[§, whereas we will show that it is quite important for color- 
suppressed modes. As a by-product, we shall see that the troublesome infrared divergent 
problem encountered in the exclusive two-body decays does not occur in the semi-inclusive 
case. 

The consideration of semi-inclusive decays B — ► MX has two more complications than 
the corresponding two-body decay, b — > Mq. First, in the free quark approximation, the 
fact that B — > MX can be viewed as the free b quark decay b — > Mq is justified only in the 
heavy quark limit. For the finite b quark mass, it becomes necessary to consider the initial 
b quark bound state effect. This has been investigated recently in || using two different 
approaches which we will follow closely. Second, the 3-body decay B — > M + qx + q 2 with 
the quark content (bq 2 ) for the B meson could be important for color-suppressed modes as 
just mentioned before. Here one needs a hard gluon exchange between the spectator quark 
q 2 and the meson M in order to ensure that the outgoing q 2 is hard. 

In passing, we briefly recall that for experimental purposes a useful feature of these 
semi-inclusive decays, B —>■ M + X, with an energetic M, is that these events should have 
relatively low multiplicity M. 

The paper is organized as follows. In Sec. II we will outline the QCD factorization 
approach relevant for our purposes and then we proceed to apply it to the two-body decays 
of the b quark b —>■ Mq in Sec. III. In Sec. IV we study the hard spectator corrections 
and summarize the results for branching ratios and CP rate asymmetries. Sec. V gives our 
conclusions. 

II. QCD FACTORIZATION 

In this section we want to suggest that the idea of QCD factorization M can be extended 
to the case of semi- inclusive decays, B — > M + X, with rather energetic meson M, say 
Em > 2.1 GeV. Recall that it has been shown explicitly M that if the emitted meson M 2 is 
a light meson or a quarkonium in the two-body exclusive decay B — » M\M 2 with M\ being 
a recoiled meson, the transition matrix element of an operator O, namely (MiM 2 \0\B) , is 
factorizable in the heavy quark limit. Schematically one has |l[ 
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{M 1 M 2 \O i \B) = (AfiM 2 |Oi|S) 



fact 



Aqcd - 



+ / d£^^l^(£,U,7;)$ B (0$M>)$M>), (2.r 
JO 



where F BMl is a B — M\ transition form factor, $m is the light-cone distribution amplitude, 
and T 1 , T 7/ are perturbatively calculable hard scattering kernels. In the naive factorization 
approach, T 1 is independent of u as it is nothing but the meson decay constant. However, 
large momentum transfer to M 2 conveyed by hard gluon exchange implies a nontrivial convo- 
lution with the distribution amplitude $m 2 - The second hard scattering function T n , which 
describes hard spectator interactions, survives in the heavy quark limit when both Mi and 
M 2 are light or when M\ is light and M 2 is a quarkonium [|lj. The factorization formula 
( [2.1D implies that naive factorization is recovered in the mj — > 00 limit and in the absence 
of QCD corrections. Nonfactorizable corrections are calculable since only hard interactions 
between the (BMi) system and M 2 survive in the heavy quark limit. 

As an illustrative example of a case when QCDF is not applicable, let us mention the 
decay B° — > 7i°D°. QCDF does not work here because the emitted meson D° is heavy so 
that it is neither small (with size of order 1/Aqcd) nor fast and cannot be decoupled from 
the (Bit) system. This is also ascribed to the fact that the soft interaction between (Bit) 
and the c quark of the D° meson is not compensated by that between (Bit) and the light 
spectator quark of the charmed meson. 

As for the semi-inclusive decay B — > MX, a momentum cutoff imposed on the emitted 
light meson M, say pm > 2.1 GeV, is necessary in order to reduce contamination from the 
unwanted background and ensure the relevance of the two-body quark decay b — > Mq. For 
example, an excess of K(K*) production in the high momentum region, 2.1 < pk(k*) < 2.7 
GeV, will ensure that the decay B — > K(K*)X is not contaminated by the background b — > c 
transitions manifested as B — > D(D*)X — > K(K*)X' and that it is dominated by the quasi 
two-body decay b — > K(K*)q induced from the penguin process b — > sg* — > sqq and the 
tree process b — > uus. As we shall see shortly, the kinematic consideration here will restrict 
the possible forms of factorization for the matrix element (MX\0\B). By the same physical 
argument as in the exclusive case and by treating M% = X and M 2 = M, it is natural to 
expect that the factorization formula (|2.1|) can be generalized to the semi-inclusive decay: 



(MX\0\B) = (MX\0\B)^ ct 



l + 5> n < + 0( 



Aqcd - 



f 1 duT I (u)$ M (u)+ I' d^duT H (^u)^ B (O^M(u). (2.2) 
Jo Jo 



In comparing this equation to the exclusive case Eq. ( |2.1|) , a crucial simplification that 
has occurred is that the semi-inclusive case does not involve any transition form factor (s). 
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This attractive feature holds so long as the meson M does not contain the spectator quark 
in the initial B meson. Since lack of knowledge of these form factors is often a serious 
limitation in quantitative applications, this adds to the appeal of the semi-inclusive case. 
Note also that when the emitted meson M is a light meson or a quarkonium, the nonfac- 
torizable corrections to naive factorization are infrared safe in the heavy quark limit and 
hence calculable. However, by the same token as the B° — > 7i°D° decay, the above QCD 
factorization formula is also not applicable to B° — > D°(D°)X. The analog of B° — ► D + n~ 
in semi-inclusive decays is K° — > ti~X c . However, there is one difference between them, 
namely the hard spectator interaction proportional to the kernel T 11 vanishes in the decay 
R -> D + n-, while it survives in the 3-body decay B — > n cq, where q is the spectator 
quark of the B meson. This is because the spectator quark in the B and D bound states in 
B — D transition is very soft in the heavy quark limit, whereas the same quark q appearing 
in the 3-body decay has to be hard so that a hard gluon exchange between tt and q is needed. 
The factorizable hadronic matrix element (MX\0\B) has the general expression: 



(MX\0\B) tact = (Mlh^iXlhlB) + (X\j[\0)(M\f 2 \B), 



(2.3) 



where j[ and j 2 arise from the Fierz transformation of the operator O and the annihilation 
term (MX\ji\0)(0\j 2 \B) is suppressed by order h*QpT>l m b an d hence it will not be included 
in Eq. (|2.3j ). As stressed in ||, Eq. (|2.3j ) is not the only way for factorization; there 
are other ways of factorization, for example, (XiM\ji\Q) (X[\j 2 \B) with X% + X[ = X. 
For the aforementioned quasi-two-body decay, it seems quite plausible to expect that the 
configuration {XiM\ji\Q) (X[\j 2 \B) is dominated by {M\jx\Q) (X\j 2 \B) @; at least in the 
perturbative region, the production of additional X\ in the final state is suppressed by 
powers of a s 0. 

The effective Hamiltonian relevant for hadronic semi-inclusive B decays of interest has 
the form 



H eS (AB = 1) 



uq 



cMo u m + cMo u m 



+ v cb v* 



cq 



Cl (n)O c M + c 2 (/i)0 2 c (/i) 
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-V tb V t * q [ y ECiQj.)O i 0i) + CgMOgiji) 



\i=3 



+VcbV* g \c 1 (n)d 1 {ji) + c 2 (n)d 2 (fi) \ \ + h.c, 



(2.4) 



where q = d,s and 



Oi = {ub) v _ A {qu) v _ A , 
O c x = (cb) v _ A (qc) v _ A , 
Oi = (cb) v _ A (qu) v _ A , 



QU 

d 2 



aJv-Ai 



(u a b /3 ) v _ A (q (3 u, 
{c a bp) v _ A (q p c a ) v _ A , 

(c*h) v - A (qpu a ) v _ A , 
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W+a)i 

(2.5) 

q> q> 

3 3 

07(9) = ^^v-AE^t^^m^Ali 08(10) = 2^Mv-A XX'O^W^v-a), 

q' q' 

O g = ±- 2 m b qa^G^(l + l5 )b, 

with q' = u, d, s, (gig 2 ) v±A = 9i7m(1 =•= 75)92, 03~06 being the QCD penguin operators, 7 - 
O10 the electroweak penguin operators, and O g the chromomagnetic dipole operator. After 
the inclusion of vertex-type and penguin-type corrections, we obtain 

T=^lv ub Vl 



V2 



uq 



ai(qu) v _ A <g> (ub) v _ A + a 2 (uu) v _ A <g> (g%_ A 
+KfeK* g [ai(gc) v _ A (8) (c&) v _ A + a 2 (cc) v ,_ A ® (#) v _ 4 _ 

+«5 EWWa ® (g&) v _ A - 2a 6 £(gg' ) S+P ® (^6) s _, (2-6) 
3 

3 3 
+-a 9 J2e q '(tf<l')v-A ® (#)v-a + 2 fl ioE e ?'(9V)v +i( ® OTv-a 
<j' 9' ' 

where the symbol <g> stands for (MX [71 ® j 2 |.B) = (M|ji|0) (X|j 2 |5) or (X|j 1 |0)(M|j 2 | J B). 

The coefficients in Eq. (|2.6| ) evaluated in the naive dimensional regularization (NDR) 

scheme for 75 have the expressions: 

c 2 a s C F 

ai = Cl + N; + ^N- c C2F > 

c x a s C F 
do = Co H Ci r , 

c 4 a s Cp 

a3 = C3 + T c + ^N- c c * F > 

^ + t + £ % K + G ^ + G ^ + |l - Cl ( a7 g ^ + " I) 
fc 3 6 1 4 1 

+(c 4 + c 6 ) 22G(si) + -(eg + cio) 22eiG(si) - -c 9 [G{s q ) + G{s b ) +-} + c g G g \, 

i=u 1 i=u 1 6 ) 

a 5 - 5 + | + ^|c 6 (-F-12), (2.7) 
«6 = c 6 + ^ + ^ ^| - Cl (^0(^) + ^0'(« C ) " I) + <*[<?(«<,) + 0'( S6 ) + \] ~ 6c 5 
+(c 4 + ce) £G'( S< ) + ^(c 8 + c 10 )^e,G'( Sl ) - ^c 9 [G'(^) + G'(s 6 ) + i + c g G' 1, 

j=« " i=u ~~ ) 



6 



c 8 a s C F a 

a7 = C7 + T c + ^N- c C8{ - F - 12) -^ NcCe > 

c 7 a s C F a , 

a8 = C8 + T c - 6 ^N- c C7 ~ 6 ^ C " 

cio a s C F a 
ag = c 9 + — + — — cio F-— N c C e , 
N c An N c 97T 

c 9 a s C F a 
ai0 = Cl0+ N- c + ^N- c C9F -^ Ce > 

where C F = {N 2 - l)/(2N e ), s t = m 2 /m 2 , and \ q = V qb V* ql . 

In Eq. (|2.7p , the vertex correction in the NDR scheme is given by [I] 



with 



F = -121n— -I8 + /7, (2.8) 
m b 

//= / dx$ M (z) (3— —]nx-3iTr), (2.9) 
jo V 1 — x J 

where $ M (x) is the leading twist light-cone distribution amplitude (LCDA) of the light 
meson M. For the vector meson V, $ (x) is dominated by the longitudinal DA (&J(x)) as 
the contribution from the transverse LCDA ($±(x)) * s suppressed by a factor of myjm b . If 
the emitted meson is the J/tp meson, then one has to take into account the charmed quark 
mass effect ||[7|]: 

* ~.T/,i,, , f„l - 2x , „. „, x 2(1-2;) 



///^= / dar^fx) 3^^1nx-3z7T + 31n(l-^) +2 
JO 11 1 — X 



1 — zx 




x \ 9 , ^ 2 x 2 [ln(l — 2) — Z7rl 

' z 2 x In zx + 



[l-^(i-x)] 2 ; [i-z(i-a:)] 5 

1 1 \ ln(l — z) — in 

^ mzx — 



, , , , ,, . r (2-10) 

X) I — ZX J 1 — -2(1— X) 

where z = m 2 j^/m 2 B and r = (/J^m c )/(/j/^), with /J^ being the tensor decay con- 
stant of the defined by 

(J/>(P, A)|or^c|0) = -ifJ/4e*^P v - 4 (A) P^). (2.11) 

Note that the third line in Eq. (|2~ToD arises from the transverse polarization component of 

II 

, / f \ ill t v/n 1 in ■ - 1 1 If i. 



the J/V- Since the asymptotic form of the distribution amplitudes (a?) and is 
the same, namely 6x(l — x), we will thus assume (x) in general. It should be 



stressed that the strong phase in // or // comes from final-state rescattering due to the 
hard gluon exchange between the outgoing M and q. 

There are QCD penguin-type diagrams induced by the 4-quark operators Oj for i = 
1,3,4,6,8,9,10. The corrections are described by the penguin-loop functions G(s) and 
G'(s) = G p (s) + G a (s) given by 
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G(s) = -^ln— + 4 f dx$ M {x) f duu(l - u) ln[s - xu(l - u)\, 
G p (s) = — In — + 3 / dx $p(x) / duu(l — u) hx[s — xu(l — u)], 

TTT/f, J «/ 

G a (s\ = --\n-^ + - C dx^-^- f 1 duu(l-u) 



3 mi, 3io x Jo 

x { In a - m l - « - 1 n ' A , 2.12 

L 2 s - xm(1 — w) j 

where and are the twist-3 LCDAs of a pseudoscalar meson P, to which we will come 
back shortly. It should be stressed that the penguin-loop contribution proportional to G'(s) 
is available only if the emitted meson is of the pseudoscalar type. In Eq. ( |2.7| ) we have also 
included the leading electroweak penguin-type diagrams induced by the operators 0\ and 

<?.= (£<?(..) + £<? W - |) (2 , 3 ) 
The dipole operator 9 will give a tree-level contribution proportional to 

Gs = _ 2 f^i^M, 

Jo x 

ff. = -f£^<,)-i (2,4) 

It is well known |J that strong phases can be perturbatively generated from the penguin 
loop functions G(x) and G'(x). The virtual gluon's momentum squared k 2 = (1 — x)m 2 B is 
no longer an arbitrary parameter; it is convoluted with the emitted meson wave function 
$ M (x). 

The twist-3 DAs 4>p and <f>^ are defined in the pseudoscalar and tensor matrix elements 



p 



(P(p)|gi(0)<v T5 g 2 (x)|0> = ~!pH x 



1 - 



nil + m 2 x 2 



1 

vqp-xiPl 



Kip***-?,,**) df}e^ x <^(f}), (2.15) 







where = nip/ {mi + rn.2) is proportional to the chiral condensate. In the present paper 
we will take the asymptotic forms for twist-2 and twist-3 distribution amplitudes: 

= 6z(l - x), $J(x) = l, (x) = 6x(l - x). (2.16) 
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Several remarks are in order, (i) The coefficients a 8a and aio a [|ll]] induced by the elec- 
troweak penguin operators 08,9,io are absorbed in our case into a 6 and 04 in Eq. ( |2.7| ). The 
contributions of C e and C' e to the electroweak coefficients 07-10, which have been neglected 
in most recent literature, are numerically more important than a 8a and a Wa owing to the 
large Wilson coefficients C\ and c 2 . (ii) The scale and 75-scheme dependence of the Wilson 
coefficients is canceled by the perturbative radiative corrections. In particular it is 

easily seen that the scale dependence of a, is compensated by the logarithmic \i dependence 
in F. However, note that the coefficients and ag are scale dependent. This is because 
the hadronic matrix element of (S — P)(S + P) operator is proportional to fi^/m b which 
is also scale dependent owing to the running quark masses, the scale dependence of and 
a 8 is canceled by the corresponding one in /z x (/i)/m&(/z). We have included penguin-type 
corrections to 04,6,8- Moreover, we have included the new contributions from the twist-3 DA 
Q a (x). (iii) In the generalized factorization approach for nonleptonic decays, the nonfactor- 
ized effects are sometimes parametrized in terms of the effective number of colors N° s so 
that 

1 1 

a-u - c 2i + ,, refn c 2 i-i, a 2 i-i - c 2 i-i + . ff . c 2i . (2.17) 

Furthermore, it is assumed that (N° s )i is process independent. In the improved QCD fac- 
torization approach, Oj (see Table II) and hence (N° s )i are in general complex and they are 
process and % dependent. 



III. TWO-BODY DECAYS OF THE b QUARK 



The general decay amplitudes for b — > Pq and b — > Vq read 

A(b -> Pq) = l ^{(A t V ub V: ql - JW a VZ)fp&qy?Q- ~ 1^ 
-BV tb V*,f P q(l- l5 )b}, 



A(b 


-> Pc) 


= & A'V^fpfp c 7m (1 - 75)6, 




A(b 


-Vq) 


= ^(A'VrtVb - A p V tb V; q ,)f v m v e^q lfl (l - 


- 75)6, 


A(b 


-> Vc) 


= ^ A'V^fvmve^ q llx (l - 75)6, 




A(b^ 




= ^(IV^ - AW^VQfj^mj/^q^ 


1 - 75)6, 



where q' = d,s is not necessarily the same as q, and the superscripts t and p denote tree and 
penguin contributions, respectively. The coefficients A and B relevant for some two-body 
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hadronic b decay modes of interest are summarized in Table I. Owing to the complications 
for the r) and production, their cases will be discussed separately below. Note that the 
coefficient B proportional to fi x /m b is formally power suppressed in the heavy quark limit, 
but numerically it is important since n x /m b ~ 12AQ CD /m 6 [see a discussion after Eq. ( |3.15| )]. 
Therefore, we will keep this calculable power correction. 

The decay amplitudes of b — > rf^s, rj^d have the expressions 



A(b -> i] n s 



V ub V: s a 2 f; (l) + V cb V; s a 2 f; C) - V tb V*((a 3 -a 5 -a 7 + ajffa 



+ [a 3 + a 4 - a 5 + -(a 7 - a 9 - a 10 )]/*o 
+ (2a 3 - 2a 5 - ^a 7 + -a 9 )/^ (0 



+ (2a 6 - 



2m s {m b — m s 



■(/;c)-/;',mi-7 5 )6 



(3.2) 



and 



A(b -> r/'^) = i 



7^' 



VubV^f^ + V cb V c * d a 2 f^ - V tb V t * d ({a 3 -a 5 -a 7 + a 9 )f^ } 



+ [2a 3 + a 4 - 2a 5 + -(-a 7 + a 9 - a w )]f%, 



+ [a 3 - a 5 + -(a 7 - a 9 )]/* ( / 



^0^(1 -75)6 



+(2a 6 - a 8 ) 



2m. 



^(/,v,-r„)v»«-^ 



(3.3) 



where the decay constants of the 77 and rj' are defined by (0| 57^755 = if q ^Pfi- The 77' 
decay constants follow a two-angle mixing pattern |T 



(3.4) 



with fs and /o being the decay constants of the SU(3) octet and singlet rjs and 770, respectively: 



(0\A°\ m ) = if p„ (0\Al\r, 8 ) =ifi 



(3.5) 



Likewise, for the r/ meson 



n = — cos 6* 8 - -y=sint) , 



(3.6) 



It must be emphasized that the two-mixing angle formalism proposed in [[L^IH] applies to 
the decay constants of the 7/ and 77 rather than to their wave functions. It is found in ||13 
that phenomenologically 
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# 8 = -21.2°, 9 = -9.2°, /„//„ = 1.26, f /f n = 1.17. 
Numerically, we obtain 

L^ = 78MeV, f*- 



-112 MeV, 



rU,d 
Iff 



63 MeV, f s v , = 137 MeV. 



(3.7) 



(3i 



The decay constant /f>, defined by (0|c7 M 7 5 c|77') = if^q^, is related to the intrinsic charm 
content of the 77' and it has been estimated from theoretical calculations JT4I] and from the 
phenomenological analysis of the data on J ftp — ■> ?7 C 7, J/-0 — > 7/7 and of the 777 and 7/7 
transition form factors |S,13,15 ; it is expected to lie in the range -2.0 MeV < f°, < -18.4 
MeV. In this paper we use the values 



fc, = -(6.3 ± 0.6) MeV, /° = -(2.4 ± 0.2) MeV, 



(3.9) 



as obtained from a phenomenological analysis performed in [13 



Care must be taken to consider the pseudoscalar matrix element for rf^ — > vacuum 
transition. The anomaly effects must be included in order to ensure a correct chiral behavior 
for the pseudoscalar matrix element [16]. The results are [T7 ; j| 

™2 



</Wlo> = -^ (fy -#>). 

(/)|«75U|0) = (r]^\d l5 d\0) = r (o (v {,) \si 5 s\0), 



(3.10) 



with 16 



Trf 



Vi-fi costf + ^sinfl 



2 /| - /2 cos 6» - v^sin 6> ; 



r„ = — 



1 Pfo- fi cos ^-^2 sin^ 



2 v^^ co ^ + 72 sin ^ 

where # is the 77 — 77' mixing angle: 

7/ = 77s sin 6 + 770 cos 0, 77 = r/ 8 cos 9 — r] sin 0. 



(3.11) 



(3.12) 



We will use # 



■15.4° as determined in 13 



To proceed with numerical calculations, we need to specify some input parameters. For 
Wilson coefficients, we use the next-to-leading ones in the NDR scheme 

ci = 1.082, c 2 = -0.185, c 3 = 0.014, c 4 = -0.035, c 5 = 0.009, c 6 = -0.041, 

c 7 / a = -0.002, c 8 /a = 0.054, c 9 /a = -1.292, c w /a = 0.263, c g = -0.143, (3.13) 

with a being the electromagnetic fine-structure coupling constant. These values taken from 



Table XXII of [pl| are evaluated at \i = m^m^) 
decay constants other than 77 and 77' we employ 



4.40 GeV and = 225 MeV. For the 

MS 
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f n = 132 MeV, f K = 160 MeV, f p = 216 MeV, f K . = 221 MeV, 

f u = 195 MeV, U = 237 MeV, / J/v , = 405 MeV. (3.14) 

For the CKM matrix elements, we take \V c b\ = 0.0395 and |Kt&/Kb| = 0.085. For the 
unitarity angle 7 we use 7 = 60° as extracted from recent global CKM fits [B5|. The other 
two unitarity angles a and (5 are then fixed. 

The hadronic matrix element of the (S — P)(S + P) operator involves the quantity ^/mj 



p 9 

= £. (3.15) 

m b (n) m 6 (/x)[mi(/x) + rn 2 (n)Y 

where m q is the running quark mass of the quark q, and hence it is formally Aqcd/w& 
power suppressed in the heavy quark limit. However, numerically it is chirally enhanced. At 
jj, = m b scale, we have (see |20| for the running quark mass ratios and evolution) 



m s (m b ) = 90 MeV, m d {m b ) = 4.6 MeV, m u (m b ) = 2.3 MeV, (3.16) 

corresponding to m s (l GeV) = 140 MeV or m s (2 GeV) = 101 MeV. We will neglect the small 
flavor-SU(3) breaking for the chiral condensate and use the averaged value /z x (m b ) = 2.7 
GeV .[] As a consequence, fi x /fn b = 12AQ CD /m b yields a large chiral enhancement. Apart from 
the current quark masses appearing in the use of equations of motion, we will utilize the pole 
masses for the heavy quarks: m b = 4.8 GeV and m c = 1.5 GeV. To compute the branching 
ratios, we use the mean lifetime of the admixture of bottom particles: r = 1.564 x 10~ 12 s 
21 . 

Since the coefficients 04.6-10 involve the CKM matrix elements X UyC /Xt [see Eq. (|2.7| )1, the 
results of the coefficients tij are exhibited in Table II for b — > P q and q' = d, s. It is evident 
that 03,4,5,6 are enhanced substantially compared to the naive factorization values and have 
large imaginary parts and that a 2 in QCD factorization becomes very small. In particular, 
the calculated coefficient a2(Pq) = 0.02 — 0.08z in QCD factorization is dramatically different 
from the value 0.18 obtained in naive factorization. The smallness of imposes a serious 
problem. For example, the predicted branching ratio of b — > J ftps is too small compared 
to the experimental value B(B -> J/if)X a ) = (8.0 ± 0.8) x 10~ 3 |22|. We will come to this 
point in Sec. IV. Note that a 9 is the dominant electroweak penguin coefficient. Owing to the 
large cancellation between 03 and 05, the decay b — > <pd is electroweak penguin dominated. 
Likewise, b — ► tt°s and p°s are also dominated by electroweak penguin diagrams. 



*A value as small as fi x = 1.4 GeV is sometimes chosen in the literature. However, we favor the 
higher value given above as one can then readily account for the observed B — * Kir rates, which 
are difficult to explain in terms of the smaller fi x . 
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We are now ready to compute the decay rates using, 



r(b 


^Pq) 


T(b 


-> Pc) 


T(b 


^Vg) 


r(6 


-> Vc) 


r(6^ 


JH>q) 



Glflml f^ty^ _ jpvaW + \B V tb V*,\ 2 ) (l - 



16n q \ ml m A h j 

^§^\A^V ub V: q , - AW tb V^ f 1 + 4 - 2< 
167T y y y mf ml 

~l&r lAV * V ^ [ 1 ~^n^ + mi 



1 'I -lA^ - AW;f (l + ^ - 2^ ) . (3.17) 
1d7t y q \ mi m 4 ' 



GfIj/^p 



The expression of T(6 — > 77 ^g) is similar to that of T(b — > Pg). We will also study the 
CP-violating partial-rate asymmetry (PRA) defined by 

_ r(6 -» Mgj - rg -» Mg) 

r(6 -»• Mg) + r(6 -»• Mg) ' 1 ' ' 

The CP-averaged branching ratios and direct CP- violating partial rate asymmetries for 
some two-body hadronic b decays of interest are summarized in Table III. Compared to 
the predictions of branching ratios based on naive factorization there are three major 
modifications: (i) Decay modes ir~u, K°d, K*°d and K~u are significantly enhanced owing 
to the large penguin coefficients a 6 and a 4 . (ii) The modes 7T°cf, p°d, u>d, J ftps, J/ipd with 
neutral emitted mesons are suppressed relative to the naive factorization ones due to the 
smallness of a^- (iii) The <pd mode has a smaller rate due to the large cancellation between 
03 and 05. That is, while <pd is QCD-penguin dominated in naive factorization, it becomes 
electroweak-penguin dominated in QCD factorization. 

For the prompt 77' production in semi-inclusive decays, we find the four-quark operator 
contributions to b — > rfs can only account for about 10% of the measured result 



B(B -> r/X.) = (6.2 ± 1.6 ± 1.3±?;£(bkg)) x 10~ 4 for 2.0 < py < 2.7 GeV/c, (3.19) 

where X s is the final state containing a strange quark. One important reason is that there 
is an anomaly effect in the matrix element (^'|s75s|0) manifested by the decay constant f™, 
[see Eq. (|3 . 10|) ] . Since f™, ~ |ft [cf. Eq. (|3.8|) ], it is obvious that the decay rate of b — > rj's 
induced by the (S — P)(S + P) penguin interaction is suppressed by the QCD anomaly effect. 
If there were no QCD anomaly, one would have B(b — ► rj's) = 2.2 x 10~ 4 from four-quark 
operator contributions which are about one third of the experimental value. 

It is useful to explicitly take into account the constraints from the CPT theorem when 
computing PRA's for inclusive decays at the quark level [^J (for a review, see |2li|]). Take 
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b — > duu as an example and note that the penguin amplitude, say XtdA, can be written as 
— (X u a% + \ c al), where X q = V q bV* ql , a\ = —CiG(s u ) + • • • and a\ = —ciG(s c ) + • • • with the 
ellipses being the common terms given in Eq. ( |2.7| ) for and G(s q ) the penguin function 
with the internal quark q [see Eq. ( |2.12 )]. In general, has absorptive contributions from 
all u, d, s, c quark loops. It is clear that PRA's arise from the interference between the tree 
amplitude X u a± and the penguin amplitude X c a% and the one between X u a1 and X c a%. The 
CPT theorem implies that the "diagonal" strong penguin phases coming from the diagonal 
process q + q — > q + q will not contribute to the rate asymmetry p6| . For example, at order 
a s the interference between X u ai and X c al with the absorptive cut from the u quark loop in 
the penguin diagram does not contribute to PRA in b — > s(d)uu decays. It is easily seen that 
the compensating process for this interference is itself. Likewise, the PRA at order a 2 s due 
to the interference between two different penguin amplitudes with the uu absorptive cut in 
the penguin loop will be compensated by the one between the tree amplitude and the higher 
order penguin diagram that contains an absorptive part from the uu quark loop inset in the 
gluon propagator. Therefore, one has to disregard the penguin phases coming from G(s u ) 
and G'(s u ) for the PRA's in the decay b — > duu. Similarly, the phase of the s-loop penguin 
diagram should be dropped in rate-asymmetry calculations of b — > dss and b — > sss in order 
to be consistent with the requirement of the CPT theorem. By the same token, the strong 
"diagonal" phase in coefficients a« due to final-state hard gluon exchange [see Eq. Q2.9Q 1 will 
not contribute to rate asymmetries in quark level processes. 

The implication of the CPT theorem for PRA's at the hadron level in exclusive or semi- 
inclusive reactions is more complicated ||27|| . Consider the above example b — > duu again. 
The corresponding semi-inclusive decays of the b quark can be manifested as b — > (tt~ , p~)u 
and b — > (n°, p°, u)d at the two-body level and (7r~7r°, K°K~)u, (7r + 7r~, 7r°7r°, K + K~)d at the 
three-body level and etc. The CPT theorem no longer constrains the absorptive cut from the 
w-loop penguin diagram not to contribute separately to each aforementioned semi-inclusive 
b decay, though the cancellation between uu and cc quarks will occur when all semi-inclusive 
modes are summed over. In view of this observation, we shall keep all the strong phases in 
the calculation of direct CP violation in the individual semi- inclusive decay. 

The presence of the strong phase in the hard kernel // [Eq. (|2.9|) 1 in QCD factorization 
has several noticeable effects: (i) The rate asymmetries in the decays b — > <pd, (n ,p ,uj)s 
vanish in naive factorization because the coefficients a-i and 03,5,7,9,10 there are real. In QCD 
factorization, the large phase of a<i (see Table II) will yield large PRA's for these modes, 
(ii) The color-suppressed tree-dominated decays b — > (rr , p°, u)d have large PRA's due to 
the large imaginary part of a<i and the smallness of |a2[- To see this, we consider b — > p°d 
as an illustration. The partial rate difference Ar(6 — ► p°d) = T(b — > p°d) — T(b — ► p°d) is 
proportional to (see Table I for the amplitude) 
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AT(b - p°d) ex ImCKftV^^Kd) Im 



3 3 1 

«2(-04 + 2 fl 7 + 2 fl 9 + 2° 10 '' 



(3.20) 



Since a 2 is dominated by the imaginary part, it follows that 



/ 3 3 1 \ 

AT(b — > p°d) oc sin a lma 2 Re ( — a 4 + -a 7 + -a 9 + -a 10 J . (3.21) 

Because | V^&V^ d | ^ l^&Krfl aric ^ l a 2| is small, it is clear that b — ► p°d has a large PRA governed 
by the vertex-induced strong phase. In contrast, the rate asymmetry in penguin-dominated 
modes is largely due to the strong penguin phase. Consider the process b — > K*~u. The 
partial rate difference is 

Ar(6 -> K*-u) oc Im^^* VXV^)lm[ ai a\ + a u 4 a c 4 } (3.22) 
~ sin 7 [Imai Real ~ ciImG(s c )Reai — CiImG(s u )Rea4 — CilvaG(s c ) Rea^], 

where we have applied the relation A t a 4 = — (X u a2+X c al) as given before. Since the imaginary 
part of ai is very small (see Table II), it is evident that CP asymmetries in the penguin- 
dominated modes are governed by the strong penguin phase. 

In the present QCD factorization approach we have considered the penguin corrections to 
04,6-10 induced by not only tree operators but also QCD and electroweak penguin operators. 
For example, the parameters 07,9,10 do receive an absorptive contribution via electroweak 
penguin-type diagrams generated by tree 4-quark operators [see Eq. (|2.7|) 1. Therefore, 
they receive more penguin phases. Another important difference between QCDF and naive 
factorization is that the gluon's virtuality k 2 is no longer arbitrary; this tends to remove 
considerable uncertainties in the estimates of the CP asymmetries. These may account for 
the general differences between the present results and [Q. Note that our predictions for 
PRA's in the decays B° -> K~{K*-)X and B~ -> K°(K*°)X are in agreement with g for 
7 = 60°. (The definition of the rate asymmetry in B has a sign opposite to ours.) 



IV. SEMI-INCLUSIVE B DECAYS 



A major advantage of studying the quasi-two-body decay of the b quark is that it does 
not involve the unknown form factors and hence the theoretical uncertainty is considerably 
reduced. In order to retain this merit for semi-inclusive B decays, it is necessary to cir- 
cumvent the second matrix element term appearing in Eq. ( |2.3|) . Fortunately, this can 
be achieved by choosing a B meson whose spectator quark content is not contained in the 
outgoing meson M. For example, the counterpart of b — > 7T°d at the meson level will be 
B° s — > 7t°X as then the B s — tt transition does not contribute. In contrast, the decay 
B~ — > tt°X or — > Tt°X will involve the unwanted B — tx form factors. As stressed in 
Sec. II, it is necessary to impose a cut, say E M > 2.1 GeV for the light emitted meson, in 
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order to reduce contamination from the unwanted background and enhance the presence of 
the two-body quark decay b — > Mq. Therefore, in the absence of the bound state effect it is 
expected that, for example, T(B° — > ix~X) « T(b — > n~u) after applying the parton-model 
approximation 

E * S / (2 X P 2E k(pu, «))Mp», *)l- (4- 1 ) 



There are two more complications for semi-inclusive -B decays. First, — > MX can be 
viewed as the two-body decay b — > Mq in the heavy quark limit. For the finite b quark mass, 
it becomes necessary to consider the initial b quark bound state effect. Second, consider the 
3-body decay B — > Mq\q\ with the quark content (652) for the B meson. One needs a hard 
gluon exchange between the spectator quark q 2 and the meson M in order to ensure that 
the outgoing q 2 is hard. For exclusive two-body decays, the nonfactorizable hard spectator 
contribution is customarily denoted as [|TJ 

G F a s C F 

71 ^ w c cJl1 - (4 - 2) 

Numerically, fjj is much larger than // for exclusive decays. For the semi-inclusive case at 
hand, it has been argued that fu is subject to a phase-space suppression since it involves 
three particles in the final states rather than the two-body one for // [|J. However, we shall 
see below that it is not the case for color-suppressed decay modes. 

A. Initial bound state effect 

The initial bound state effects on branching ratios and CP asymmetries have been studied 
recently in using two different approaches: the light-cone expansion approach and the 
heavy quark effective theory approach. We will follow [f| to employ the second approach 
which amounts to modifying the decay rates by 



t(b - px) = - A p v tb v: q ,\ 2 m + \Bv tb v; ql \% 




T(B - VX) = C lllt±\A t V ub V: ql - AFV^W [ 1 + ^ - 2^ ] , (4.3) 

c-2 



r(B -> j/^x) = GFf ^f lb \A t v d y; q - A?v tb v; q \% [ 1 + 

where 







m% 








m j/1> 






z 4 



7 /j,q 53 ^\ / ii G K , , A 
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and 



(B\hG^ u h\B) 2 (B\h(iD ± ) 2 h\B) 



with Z}^ = £> M — v^v ■ D, where v is the four- velocity of the B meson. The nonperturbative 
HQET parameter p? G is fixed from the B*—B mass splitting to be 0.36 GeV 2 . Following || we 



use p 2 , = 0.5 GeV 2 , which is consistent with QCD sum rule and lattice QCD calculations [p8 



Compared to the two-body decays b — > Mq shown in Table III, we see that the branching 
ratio of B — > PX and B — > VX owing to bound state effects is reduced by a factor of 
(5 ~ 10)% and 17%, respectively, while the CP asymmetry remains intact for VX decays 
and for most of PX modes. 



B. Nonfactorizable hard spectator interactions 

We now turn to the hard spectator interactions in the 3-body decay B(p B ) — > M(j>m) + 
Qii.Pi) + ^2(^2) with a hard gluon exchange between the spectator quark q 2 and the meson 
M. A straightforward calculation yields 

A spect (B - P gi q 2 )= ^ ^ ^ ^ /p/B«K fe KV - 

x f 1 d ^ P ^ f 1 dp^^ (pUil^-l^)q2 m B q l (l+^)q 2 \ 



h £ Jo p \ p P ■ p 2 p B - P2 

A spcct (B - V qi q 2 )= -i^ ^ ^ ^ m v f v f B (Ai p V ub V: q , - A%V tb V; ql ) 



x 



1 ,1 
d£ dp ■ 
Jo 



Jo £p(PB-P2)(pvP2 + £mlr/2) 
x |(Pb ■ P2 - iPB • p v )e*^qii^{l - 75)92 - • Pb)9i?V(1 ~ 7s)92 

-m B (e* ■ p 2 )qi(l + 75)92 - €m B qip'vf(l + 75)92 j, 
A spect (S^ J/^fc)=-i^|i^^ mj/ ^ d^J o dp 

x $J/ ^(0$f(p) I [j9b . p2 _ ^ _ r)pB . ]e .p (1 _ ^ 

£p(Pb ■P2){Pj/i, -P2) { 
-(£ - r)(e* ■ p B )qii>j/^{ 1 ~ 75)^2 - ■ £2)91(1 + 75)92 

-^-r)m B qi^j/^*(l + 75)92! , (4.6) 

where r = (m c fj^)/(mj/^fj/^). In deriving the above equation, we have applied the on- 
shell conditions 9iA = 0, ^ 2 <?2 = 0, the approximation p ps [see the discussion after Eq. 
( |1.11| )1 and the i? meson wave function 0: 
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(0\q a (x)bp(0)\B(p))\ x+=x±=0 = -^[(^ + m B ) 75 ]^^ 1 dpe-^[$f(p) + ^_<&f (p)] 7a , (4.7) 

with n_ = (1,0,0,-1). The expressions for A 1 ^ and A 1 ^ can be obtained from that of 
A t,p and A t,p (see Table I) respectively by replacing the coefficient 0,21 (a2i-i) by the Wilson 
coefficient c<n-\ [c-ii)- 111 passing, we note that, in contrast to f|, the 2-body decay b — > Mq\ 
and the 3-body decay B —>■ Mq\q2 do not interfere with each other to give contributions to 
B -> MX. 



Eq. ( |4.6|) can be applied to two-body exclusive decays. Consider B — > Ktc as an example 
and this amounts to having P = K and a pion from q\q2 . Hence P2 = fjPn, where fj is the 
momentum fraction of the antiquark in the pion. It follows from (f4.6|) and ( |2.15 ) that 



A spect (B - Kn) oc f 1 ? $f (p) /' f $^(0 ? 

Jo p io ( Jo 1] 



m b p 



which was first obtained in [ 29] . It is evident that the terms proportional to in Eq. 



give twist-3 contributions. Since the twist-3 distribution amplitude <&p(fj) ~ 1, it does not 
vanish at the endpoints. Consequently, there is a logarithmic divergence of the fj integral 
which implies that the spectator interaction in B — > decay is dominated by soft gluon 
exchange between the spectator quark and quarks that form the emitted kaon, indicating that 
QCD factorization breaks down at twist-3 order. The above-mentioned infrared divergent 
problem does not occur in the semi-inclusive decay, however. 

The decay distribution due to hard spectator interactions is given by 

1 1 r 

dT spect = -—3 ^—3- J \A spect \ 2 dm 2 12 dm 2 23 , (4.9) 

or 

dr spect _ 1 1 r . \2 A ^2 



dE M (2vr) 3 16m| 



Inspect I dm 2 23} (4.10) 



where Em is the energy of the outgoing meson M, and m?- = (p, +Pj) 2 with p 3 = pm- For a 
given Em, the range of m| 3 is fixed by kinematics. In order to enhance the possibility that 
B — > MX originates from a quasi-two-body decay, we impose the energy cutoff E M > 2.1 
GeV for light mesons and E M > 3.3 GeV for the J/ip. 



C. Results and discussions 



To proceed we apply the initial bound state effect to hard spectator interactions and use 
the B meson wave function 



$f(p) = iV B p 2 (l-p)'exp 



1 / pm B \ 2 

2 V uj b ) 



(4.11) 
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with ujb = 0.25 GeV and N B being a normalization constant. This B meson wave func- 
tion corresponds to X B = 303 MeV defined by Jq dp & B (p) / p = m B /\ B This can be 
understood since the B meson wave function is peaked at small p: It is of order ttib/Aqcd 
at p ~ Aq CD /m B . Hence, the integral over <p B (p)/p produces an m B /Aq C D term. As for 
the parameter r defined after Eq. ( J2.1U ), it is equal to 1/2 in the heavy quark limit as- 
suming fj^ = fj/^. The results of calculations are shown in Table IV. We see that the 
tree-dominated color-suppressed modes (n , p°,u)X s , (fiX, J/ipX s , J/ipX and the penguin- 
dominated mode ujX S s are dominated by the hard spectator corrections. In particular, the 
prediction BiB — > J/ipX s ) = 9.6 x 10~ 3 is in agreement with the measurement of a direct 
inclusive J/ip production: (8.0 ± 0.8) x 1(T 3 by CLEO |2f and (7.89 ± 0.10 ± 0.34) x 1(T 3 
by BaBar ||30|| . This is because the relevant spectator interaction is color allowed, whereas 
the two-body semi-inclusive decays for these modes are color-suppressed. As a consequence, 
nonfactorizable hard spectator interactions amount to giving a 2 a large enhancement. 

It is instructive to compare the enhancement of the J/tp production in exclusive and 
semi-inclusive decays. The hard spectator diagrams denoted by fn have been included in 
the leading-twist order calculations and it is found that a 2 ( J/ipK) ~ 0.06 — 0.05i [^,0. 
Therefore, the real part of ci2(J/il>K) is enhanced by which is numerically much larger 
than //, but it is still too small compared to the experimental value 0.26 ± 0.02 ||31|| . It has 
been shown recently that to the twist-3 level, the coefficient a 2 {J/ipK) is largely enhanced 
by the nonfactorizable spectator interactions arising from the twist-3 kaon LCDA (p^ , which 
are formally power-suppressed but chirally, logarithmically and kinematically enhanced ||. 
The major theoretical uncertainty is that the infrared divergent contributions there should 
be treated in a phenomenological way. In this work we found that it is the same spectator 
mechanism responsible for the enhancement observed in semi-inclusive decay B — > J/ipX s , 
and yet we do not encounter the same infrared problem as occurred in the exclusive case, 
and terms proportional to in Eq. (|4.6|) are not power suppressed, rendering the present 
prediction more reliable and trustworthy. It is conceivable that infrared divergences residing 
in exclusive decays will be washed out when all possible exclusive modes are summed over. 

It is also interesting to notice that after including the spectator corrections, the branch- 
ing ratios and PRA's for the color-suppressed modes B^ s — > (ir°, p°, uj)X s , B~ — > <f)X are 
numerically close to that predicted in j|] based on naive factorization (see Table III). Note 
that the large CP asymmetries in b — > (tt°, p°, u)d decays (see Table III) are washed out to 
a large extent at hadron level by spectator interactions. By contrast, the nonfactorizable 
spectator interaction is in general negligible for penguin dominated (except for ujX sS ) or 
color-allowed tree dominated decay modes. The channels (B~ , 5°) — > (tt°, p°, u, 4>)X are not 
listed in Table IV as they involve the unwanted form factors. For example, B~ — > n°X con- 
tains a term a2F Bn and B° — > tt°X has a contribution like a^F Bn . Hence, the prediction of 
(B~ , B°) — > 7r°X is not as clean as B~° s — > 7r°X 5 . Nevertheless, the former is also dominated 
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by spectator interactions and is expected to have the same order of magnitude for branching 
ratios as the latter. 

Owing to the presence of B — 77(77') form factors, the decays B — > (77, r)')X are also 
not listed in Table IV. However, we find that the hard spectator corrections to the prompt 
i]' production in semi-inclusive decays are very small and hence the four-quark operator 
contributions to b — > rj's can only account for about 10% of the measured result, Eq. (|3.19f) . 
Evidently this implies that one needs a new mechanism (but not necessarily new physics) 
specific to the rj'. It has been advocated that the anomalous coupling of two gluons and 77' 
in the transitions b —>■ sg* followed by g* — ► rj'g and b — > sg*g* followed by g*g* — > 77' may 
explain the excess of the 77' production |32||l7j . An issue in this study is about the form-factor 
suppression in the 77' — g* — g* vertex and this has been studied recently in the perturbative 
QCD hard scattering approach |33|]. At the exclusive level, it is well known that the decays 
B ± — > r)'K ± and 5° — > r)'~K° have abnormally large branching ratios |HJ] . In spite of many 
theoretical uncertainties, it is safe to say that the four-quark operator contribution accounts 
for at most half of the experimental value and the new mechanism responsible for the prolific 
77' production in semi-inclusive decay could also play an essential role in B — > rj'K decay. 

From Table IV it is clear that the semi-inclusive decay modes: B° s — > (tt°, p°, u)X s , p°X S s, 
B -> (K~X, K*-X) and Br -> (K°X S , K*°X S ) are the most promising ones in searching for 
direct CP violation; they have branching ratios of order 10~ 6 — 10~ 4 and CP rate asymmetries 
of order (10 — 40)%. Note that as shown in Eqs. ( |3.21| ) and ( 3.22 ), a measurement of 
partial rate difference of B° s — > (7r°, p°, u)X s and B~ — > (K°X S , K*°X S ) will provide useful 
information on the unitarity angle a, while B s — > p°X sS and B — > (K~X, K*~X) on the 
angle 7. To have a rough estimate of the detectability of CP asymmetry, it is useful to 
calculate the number of B — B pairs needed to establish a signal for PRA to the level of 
three statistical standard deviations given by 

< = A^< (4 ' 12 > 

where A is the PRA, Br is the branching ratio and e e g is the product of all of the efficiencies 
responsible for this signal. With about 1 x 10 7 BB pairs, the asymmetry in K*~ channel 
starts to become accessible; and with about 7 x 10 7 BB events, the PRA's in the other 
modes mentioned above will become feasible. Here we assumed, for definiteness, e c g = 1 and 
a statistical significance of 3a as in Eq. ( 4.12|) . Currently BaBar has collected 23 million 
BB events, BELLE 11 million pairs and CLEO 9.6 million pairs. It is conceivable that 
CP asymmetries in semi-inclusive B decays will begin to be accessible at these facilities. 
Likewise, PRA's in semi-inclusive B s decays may be measurable in the near future at the 
Fermilab's Tevatron. 

It is interesting to note that the decays B° s — > (tt°, p ,u)X S g and B~ — > <pX are 
electroweak-penguin dominated. Except for the last channel, they have sizable branching 
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ratios and two of them have observable CP asymmetries. A measurement of these reactions 
will provide a good probe of electroweak penguins. 

Finally, it is useful to discuss briefly the theoretical uncertainties one may have in the 
present approach for semi-inclusive decays: the b quark mass, the annihilation diagram and 
the distribution amplitude of the meson. We have assumed a pole mass for the b quark 
to compute the decay rates of b — > Mq and B — > MX which are proportional to m\. In 
principle, this uncertainty in can be reduced by normalizing the semi-inclusive hadronic 
rate to the semi-leptonic one. Since the latter is of 5th power in m&, the uncertainty is only 
slightly alleviated. The annihilation topology is power suppressed in the heavy quark limit 
and is conventionally assumed to be small. However, it is conceivable that power corrections 
due to the annihilation diagrams could be important for penguin-dominated semi-inclusive 
decays such as B~ — > K°(K*°)X and B° — > K~ (K*~)X . As for the LCDAs of the meson, 
we have assumed the asymptotic form ( |2.16| ) for the leading-twist LCDA, which is suitable 
for light mesons but probably not so for the heavy meson J/ifi. Due to SU(3) symmetry 
breaking, the realistic kaon wave function should exhibit a slight asymmetry in 1 — 2x. Also 
the distribution amplitude of the B meson is not well understood; phenomenologically, the 
parameter ujb [see Eq. ( [4.111) 1 or Ab is not well fixed. 

V. CONCLUSIONS 

We have systematically investigated semi-inclusive B decays B — > MX within the frame- 
work of QCD-improved factorization. The nonfactorizable effects, such as vertex-type and 
penguin-type corrections to the two-body b decay, b — > Mq, and hard spectator corrections 
to the 3-body decay B — > Mq\q\ are calculable in the heavy quark limit. QCD factorization 
is applicable when the emitted meson is a light meson or a charmonium. 

There are two strong phases in the QCD factorization approach: one form final-state 
rescattering due to hard gluon exchange between M and X, and the other from the penguin 
diagrams. We have discussed the issue of the CPT constraint on partial rate asymmetries. 
The strong phase coming from final-state rescattering due to hard gluon exchange between 
the final states M and X [see Eq. (|2.9|) 1 can induce large rate asymmetries for tree-dominated 
color-suppressed modes (ir°, p°, uj)X s . The predicted coefficient a<i in QCD factorization 
is very small compared to naive factorization. Consequently, the color-suppressed modes 
(7T°, p°, u)X s , (fiX and J/ifiX s , J/ipX are very suppressed. Fortunately, the nonfactorizable 
hard spectator interactions in B — > Mqxq 2 , though phase-space suppressed, are extremely 
important for the aforementioned modes. In fact, they are dominated by the hard spectator 
corrections. This is because the relevant hard spectator correction is color allowed, whereas 
the two-body semi-inclusive decays for these modes are color-suppressed. Our prediction 
BiB — > J/ifiX s ) = 9.6 x 10 -3 is in agreement with experiment. Contrary to the exclusive 
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hadronic decay, the spectator quark corrections here are not subject to the infrared divergent 
problem, rendering the present prediction more clean and reliable. 

Owing to the destructive QCD anomaly effect in the matrix element of pseudoscalar 
densities for ?/-vacuum transition, the four-quark operator contribution to b — > rj's is too 
small to explain the observed fast production. It is evident that a new mechanism such as 
the anomalous coupling of two gluons with the rj' is needed in order to resolve the 7/ puzzle. 

The semi-inclusive decay modes: ~B^ S — > (tt°, p°, uj)X s , p°X sS , B° — > (K~X, K*~X) and 
B~ — > (K°X S , K*°X S ) are the most promising ones in searching for direct CP violation; they 
have branching ratios of order 10~ 6 — 10~ 4 and CP rate asymmetries of order (10—40)%. With 
about 7x 10 7 BB pairs, CP asymmetries in these modes may be measurable in the near future 
at the BaBar, BELLE, CLEO and Tevatron experiments. The decays B° s — > (it , p°, uo)X s 
and B~ — > <pX are electroweak-penguin dominated. Except for the last mode, they in general 
have sizable branching ratios and two of them have observable CP asymmetries. The above- 
mentioned reactions will provide good testing ground for the standard model and a good 
probe for electroweak penguins. 
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TABLE I. The coefficients A t (A t ), Ap(Ap) and B (in units of fJ, x /m b ) denned in Eq. (|1|) for 
some modes of interest. 



Mode 




At I A~t\ 
A. [A ) 


API AV\ 

A y [A y ) 


B(n x /m b ) 


o - 


-> 7T U 


7 
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a 1 
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b - 


-> p~u 
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Oi 


a 4 + aio 




7 

b - 


.1 


j 
a 


a 2 /V2 


(-04 - 2 a 7 + 2 a 9 + 2 a 10)/v2 


„ / AT 


7 

o - 


^ 

-» p a 


d 


/ AT 

a 2 /V2 


/ 1 3 1 3 1 1 \ / AT 

(-a 4 + §a 7 + fa 9 + 1 aiQ)/\J2 




7 

o - 


->■ a; a 


7 

a 


/ AT 

a 2 /V2 


To 1 1 1 1 ( 1 M / /o 

[2a 3 + a 4 + 2a 5 + ^{a-j + a 9 - ai )J/V2 




7 

b - 


(pa 


7 

a 




a3 + as - 5(07 + a 9 ) 




7-, 

b - 


-* 7/ d 


7 

a 




see text 




7 
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-> T] a 


j 
a 




see text 




b - 
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7 

a 
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-> p~c 


d 
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u 
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j 
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,7 
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see text 
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see text 
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a 3 + as + a 7 + a 9 
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d 


«2 
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TABLE II. Numerical values for the coefficients ctj (in units of 10 4 for 03, • • • , aio) in QCD 
factorization and in naive factorization for b — > and q' = d,s . 



QCD factorization 





g' = d 


q' = s 


naive factorization 


ai 


1.05 + 0.0H 


1.05 + 0.01* 


1.02 


a 2 


0.02 - 0.08i 


0.02 - 0.08i 


0.18 


a 3 


74 + 26i 


74 + 26i 


23 


04 


-317 - 29i 


-353 - 58i 


-303 


a 5 


-67 - 30i 


-67 - 30i 


-8 


a 6 


-400 - 27i 


-440 - 45i 


-380 


a 7 


-0.29 - 0.65i 


-0.89 - 1.13* 


0.80 


a 8 


3.43 - 0.36i 


3.22 - 0.46i 


3.89 


a 9 


-92.3 - 2.3z 
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-87.9 


aio 


0.8 + 6.6i 


0.6 + 6.4* 


-12.2 
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TABLE III. CP-averaged branching ratios and partial-rate asymmetries for some two-body 
hadronic b decays. For comparison, the predicted branching ratios and rate asymmetries (in abso- 
lute values for 7 = 60°) based on naive factorization [4] are given in parentheses. 



Mode 


BR 








PRA(%) 


b — > ir~u 


1.5 x 


10~ 4 


(1.3 x 


10~ 4 ) 


-2 (7) 
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(5.8 x 


10- 6 ) 


-97 (34) 
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TABLE IV. C-P-averaged branching ratios and partial-rate asymmetries for some semi-inclusive 
hadronic B decays with Em > 2.1 GeV for light mesons and Ejj^ > 3.3 GeV for the J/ip. 
Branching ratios due to hard spectator interactions in the 3-body decay B — ► Mq\q2 are shown in 
parentheses. Here X denotes a final state containing no (net) strange or charm particle, and X q 
the state containing the quark flavor q. 



Mode 


BR 












PRA(%) 
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